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1 Analogue Gravity: An Overview 

1.1 The concept of an effective space-time metric 

Curved space-times are familiar from Einstein's theory of gravitation [l], 
where the metric tensor g^^, , describing distances in a curved space-time with 
local Lorentz invariance, is determined by the solution of the Einstein equa- 
tions. A major problem for an experimental investigation of the (kinematical 
as well as dynamical) properties of curved space-times is that generating a 
significant curvature, equivalent to a (relatively) small curvature radius, is a 
close to impossible undertaking in manmade laboratories. For example, the 
effect of the gravitation of the whole Earth is to cause a deviation from flat 
space-time on this planet's surface of only the order of 10~^ (the ratio of 
Schwarzschild and Earth radii). The fact that proper gravitational effects are 
intrinsically small is basically due to the smallness of Newton's gravitational 
constant G = 6.67 x 10~^^ m^kg~^sec~^. Various fundamental classical and 
quantum effects in strong gravitational fields are thus inaccessible for Earth- 
based experiments. The realm of strong gravitational fields (or, equivalently, 
rapidly accelerating a reference frame to simulate gravity according to the 
equivalence principle), is therefore difficult to reach. However, Earth-based 
gravity experiments are desirable, because they have the obvious advantage 
that they can be prepared and, in particular, repeated under possibly different 
conditions at will. 

A possible way out of this dilemma, at least inasmuch the kinematical 
properties of curved space-times are concerned, is the realization of effective 
curved space-time geometries to mimic the effects of gravity. Among the most 
suitable systems are Bose-Einstein condensates, i.e., the dilute matter- wave- 
coherent gases formed if cooled to ultralow temperatures, where the critical 
temperatures are of order Tc ^ I uK • • • 1 ^K; for reviews of the (relatively) 
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recent status of this rapidly developing field see IH[31[3]. In what follows, it 
will be of some importance that Bose-Einstein condensates belong to a special 
class of quantum perfect fluids, so-called superfluids [5]. 

The curved space-times we have in mind in the following are experienced 
by sound waves propagating on the background of a spatially and temporally 
inhomogcneous perfect fluid. Of primary importance is, first of all, to real- 
ize that the identification of sound waves propagating on an inhomogcneous 
background, which is itself determined by a solution of Euler and continuity 
equations, and photons propagating in a curved space-time, which is deter- 
mined by a solution of the Einstein equations, is of a kinematical nature. That 
is, the space-time metric is fixed externally by a background field obeying the 
laws of hydrodynamics (which is prepared by the experimentalist), and not 
self-consistently by a solution of the Einstein equations. 

As a first introductory step to understand the nature of the kinematical 
identity, consider the wave equation for the velocity potential of the sound 
field (j), which in a homogeneous medium at rest reads 



where Cg is the sound speed, which is a constant in space and time for such 
a medium at rest. This equation has Lorentz invariance, that is, if we replace 
the speed of light by the speed of sound, it retains the form shown above 
in the new space-time coordinates, obtained after Lorentz-transforming to a 
frame moving at a constant speed less than the sound speed. Just as the 
light field in vacuo is a proper relativistic field, sound is a "relativistic" field0 
The Lorentz invariance can be made more manifest by writing equation ([T]) 
in the form Dcf) = ri'^'^di^i,d,j(f> = 0, where 77^" = diag(l, — 1, — 1, — 1) is the 
(contravariant) flat space-time metric (we choose throughout the signature 
of the metric as specified here) , determining the fundamental light-cone- like 
structure of Minkowski space [5] ; we employ the summation convention over 
equal greek indices ^,v, - ■ ■ . Assuming, then, the sound speed = Cs{x,t) 
to be local in space and time, and employing the curved space-time version 
of the 3+ ID Laplacian □ [1], one can write down the sound wave equation in 
a spatially and temporally inhomogcneous medium in the generally covariant 



Here, g = dct[g^^] is the determinant of the (covariant) metric tensor. It is to 
be emphasized at this point that, because the space and time derivatives are 

^ More properly, we should term this form of Lorentz invariance pseudorelativistic 
invariance. We will however use for simplicity "relativistic" as a generic term if 
no confusion can arise therefrom. 



form [HE] 




(2) 
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covariantly transforming objects in the primary object in the condensed- 
matter identification of space-time metrics via the wave equation Q is the 
contravariant metric tensor g'"' [Hj . In the condensed-matter understanding of 
analogue gravity, the quantities g^^ are material- dependent coefficients. They 
occur in a dispersion relation of the form g'^^k^k^, = 0, where fc^ = {uj/cs,k) 
is the covariant wave vector, with hk the ordinary spatial momentum (or 
quasi- momentum in a crystal). 

The contravariant tensor components g^'' for a perfect, irrotational liquid 
turn out to be [TUHlfTU] 

(3) 



where 1 is the unit matrix and Ac a space and time dependent function, to 
be determined from the equations of motion for the sound field (see below). 
Inverting this expression according to g'^'^g^ = ^^a, to obtain the covariant 
metric g^u, the fundamental tensor of distance reads 

where the line element is ds^ = g^^dx^dx'^ . This form of the metric has been 
derived by Unruh for an irrotational perfect fluid described by Euler and con- 
tinuity equations [7]; its properties were later on explored in more detail in 
particular by M. Visser [8]. We also mention that an earlier derivation of Un- 
ruh's form of the metric exists, from a somewhat different perspective; it was 
performed by Trautman [TU]. The metric belongs to the Painleve-Gullstrand 
class of metrics, historically introduced in Refs. pT|. 

The conformal factor Ac in @ depends on the spatial dimension of the 
fluid. It may be unambiguously determined by considering the effective action 
of the velocity potential fluctuations above an inhomogeneous background, 
identifying this action with the action of a minimally coupled scalar held in 
D + 1-dimensional space-time 



2g 

d^+ixy^g^''a,,<^a,0, (5) 



where the prefactor 1/g in front of the square brackets in the first line is 
identified with the compressibility of the (barotropic) fluid, 1/g = dihi p) / dp, 
with p the pressure and p the mass density of the fluid; we assume here and 
in what follows that g is a. constant independent of space and time so that 
5Pi a-s valid for a dilute Bose gas (see the following subsection). Using 
the above identiflcation, it may easily be shown that the conformal factor 
is given by Ac = (cs/g)^^^^"^"* = (p/ff)^^*-^"^-*), while the square root of the 
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negative determinant is y/—g = Cs{cs/ g)^^'^^^^^^^ — {p^ / gY^''^^'^^ ■ The case 
of one spatial dimension (_D = 1) is special, in that the conformal invariance 
in two space-time dimensions implies that the classical equations of motion 
are invariant (take the same form) for any space and time dependent choice 
of the conformal factor A^, explaining the singular character of the conformal 
factor at the special value D = 1. 

The line element ds^ = g^^dx^dx'^ gives us the distances travelled by the 
phonons in an effective space-time world in which the scalar field "lives". 
In particular, quasiclassical (large momentum) phonons will follow /i^/ii-like, 
that is, here, sownrf-like geodesies in that space-time, according to ds^ = 0. 
Noteworthy is the simple fact that the constant time slices obtained by setting 

= in the line element are conformally flat, i.e. the quasiparticle world 
looks on constant time slices like the ordinary (Newtonian) lab space, with 
a Euclidean metric in the case of Cartesian spatial co-ordinates we display. 
All the intrinsic curvature of the effective quasiparticle space-time is therefore 
encoded in the metric tensor elements goo and goi [12] . 



1.2 The metric in Bose-Einstein condensates 



We assumed in Eq. ([5]) that the compressibility l/g is a constant. This entails 
that the (barotropic) equation of state reads p = \gp^. We then have, in the 
microscopic terms of the interaction between the particles (atoms) constitut- 
ing the fluid, a contact interaction (pseudo-)potential, V{x — x') = g5{x — x'). 
This is indeed the case for the dilute atomic gases forming a Bose-Einstein 
condensate. Well below the transition temperature, they are described, to 
lowest order in the gas parameter (pa^)^/^ ^ 1 [where as is the s-wave scat- 
tering length; assumed positive] by the Gross-Pitaevskii mean-fleld equation 
for the order parameter ^ = (ft) , the expectation valu^ of the quantum field 
operator W [3]: 



d 

ih—9{x,t) 



2 



H 

— A + Vtr^p{x,t) + g\1'{x,t)\' 
zm 



Wix.t). (6) 



Here, Vtrap denotes the one-particle trapping potential and the coupling con- 
stant g is related to the s-wave scattering length via g = ATTh?as/m (in three 
spatial dimensions). The Madelung transformation decomposing the complex 
field into modulus and phase reads 'F = ^exp[j0], where p yields the con- 
densate density and (f> is the velocity potential. It allows for an interpretation 
of quantum theory in terms of hydrodynamics [13j . Namely, identifying real 



^ Observe that {4') ^ breaks particle number conservation (the global U(l) in- 
variance) . We will come back to this point in section [2] where we introduce a 
particle-number-conserving mean-field ansatz for the full quantum field operator, 
which has the number-conserving property that {>!'} = 0, and for which the order 
parameter therefore does not equal {'I'). 
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and imaginary parts on left- and right-hand sides of ^ , respectively, gives us 
the two equations 



The first of these equations is the Josephson equation for the superfluid phase. 
This Josephson equation corresponds to the Bernoulli equation of classical 
hydrodynamics, where the usual velocity potential of irrotational hydrody- 
namics equals the superfluid phase times h/m, such that v = hV(j)/m. 
The latter equation implies that the flow is irrotational save on singular lines, 
around which the wave function phase (j) is defined only modulo 27r. Therefore, 
circulation is quantized [Hj, and these singular lines are the center lines of 
quantized vortices. The usual classical terms in Eq. ([7]) are augmented by the 
"quantum pressure" pg = —-^(Ay/p)/^. The second equation ([8]) is the 
continuity equation for conservation of particle number, i.e., atom number in 
the superfluid gas. The dynamics of the weakly interacting, dilute ensemble 
of atoms is thus that of a perfect Euler fluid with quantized circulation of sin- 
gular vortex lines. This is true except for regions in which the density rapidly 
varies and the quantum pressure term »q becomes relevant, which happens 
on scales of order the coherence lengtflfl ^ — h/y/gmp where p is a constant 
(asymptotic) density far away from the density-depleted (or possibly density- 
enhanced) region. This is the case in the density-depleted cores of quantized 
vortices, or at the low-density boundaries of the system. The quantum pres- 
sure is negligible outside these domains of rapidly varying and/or low density. 
The whole armoury of space-time metric description of excitations, explained 
in the last section, and based on the Euler and continuity equations, is then 
valid for phonon excitations of a Bose-Einstein condensate, with the space- 
time metric as long as we are outside the core of quantized vortices and 
far from the boundaries of the system, where both the flow is irrotational and 
the quantum pressure is negligible. 

We mention here in passing that the form ^ of the wave equation is 
valid in quite general physical contexts. That is, a generally covariant curved 
space-time wave equation can be formulated not just for the velocity pertur- 
bation potential in an irrotational Euler fluid, for which we have introduced 
the effective metric concept. If the spectrum of excitations (in the local rest 
frame) is linear, w = Cprop^, where Cpj-op is the propagation speed of some 
collective excitation, the statement that an effective space-time metric exists 
is true, provided we only consider wave perturbations of a single scalar field 
constituting a fixed classical background: More precisely, given the generic 
requirement that the action density £ is a functional of <P and its space-time 

^ Note that the coherence or "healing" length is also frequently defined in the 
literature, see e.g. [1], with an additional factor of l/\/2, i.e., as = h/y/2gmp. 



-p + V • (p-u) = 0. 



2m yfp 



(7) 



(8) 
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derivatives 9^?^, i.e. £ ~ C[<P, d^(f\, the fluctuations (f) = d^P around some clas- 
sical background solution <Po of the Euler-Lagrange equations always satisfy 
a wave equation of the form of Eq. with a possible additional scalar po- 
tential term |15| comprising, for example, a mass term for the scalar field. As 
a consequence, the effective metric description also applies, inter alia, to the 
quasiparticle excitations around the gap nodes in the superfluid '^le-A |16j . 
photon propagation in dielectrics [i7\ . and surface waves in shallow water 

m- 

The analogy between photons propagating on given (curved) space-time 
backgrounds and phonons in spatially and temporally inhomogeneous super- 
fluids or, more generally, quantized quasiparticles with linear quasiparticle 
dispersion in some background, allows us to apply many tools and methods 
developed for quantum fields in curved space-times [19] . We can therefore 
conclude that the associated phenomena occur (provided the fundamental 
commutation relations of these quantum fields are fulfilled [10]). Among these 
phenomena are Hawking radiation |21j , the Gibbons-Hawking effect in dc Sit- 
ter spacetime [22], and cosmological particle production [23] [21]. Furthermore, 
cosmic inflation and the freezing-in of quantum vacuum fluctuations by the 
horizon crossing of the corresponding modes may be simulated |25| . A com- 
prehensive recent review of the subject of analogue gravity in its broadest 
sense is given in [26] . 

1.3 Pseudo-Energy-Momentum Tensor 

An important quantity characterizing the dynamics of the fleld is the 
pseudo-energy-momentum tensor, cf. |27| . Since the equation of motion for 
the scalar mode, V^g^''V,y(/) = 0, where denotes the space-time-covariant 
derivative, is equivalent to covariant energy-momentum balance, expressed by 
VpT'"' _ ^j^jj classical pscudo-energy-momentum tensor reads [27] 

T^., = (a^0)(a,0) - i g,Adpci^)idA)g'^ ■ (9) 

We already stressed that the identiflcation of fleld theoretical effects in curved 
space-time by analogy (including the existence of the pseudo-energy-mo- 
mentum tensor), is of a kinematical nature. An important question concerns 
the dynamics, that is, the backreaction of the quantum fluctuations of the 
scalar fleld onto the classical background. Extending the analogy to curved 
space-times a bit further, one is tempted to apply the effective-action method 
(see, e.g., [TH] and [55]). In the effective-action method, one integrates out 
fluctuations of the quantum flelds to one-loop order, and then determines 
the expectation value of the energy-momentum-tensor by the canonical iden- 
tiflcation SAes/Sg'^" = Since the dependence of the effective 
action AcS on the degrees of freedom of the background 77 enters via the ef- 
fective metric g^'' = g^'^^if), one flnds the backreaction contribution to the 
equations of motion of the 77 by differentiation of the effective action according 
to 
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5j] (Jg^"^ Sri 2 ^ ° \ / Sr] 

The precise meaning of the expectation value of the pseudo-energy-momentum 
tensor, (T^i/), is difhcuh to grasp in general, due to the non- uniqueness of the 
vacuum state in a complicated curved space-time background and the ultra- 
violet (UV) renormalization procedure. Adopting a covariant renormalization 
scheme, the results for (T^i/) can be classified in terms of geometrical quanti- 
ties (cf. the trace anomaly [19]). 

However, in calculating the quantum backreaction using the effective- 
action method, one is implicitly making two essential assumptions: first, that 
the leading contributions to the backreaction are completely determined by 
the effective action in Eq. ([5]), and, second, that deviations from the low-energy 
effective action at high energies do not affect the (renormalized) expectation 
value of the pseudo-energy-momentum tensor, (T^^). Since the effective covari- 
ance in Eq. ^ is only a low-energy property, the applicability of a covariant 
renormalization scheme is not obvious in general. In the following, we criti- 
cally examine the question of whether the two assumptions mentioned above 
are justified, e.g., whether (T^jy) completely determines the backreaction of 
the linearized quantum fluctuations. 

The remainder of these lecture notes, which is based on the results of 
the publication [55], is organized as follows. In section [21 we give a brief 
introduction to Bose-Einstein condensates, and introduce a particle- number- 
conserving ansatz for the field operator separated into condensate, single- 
particle and multi-particle excitation parts. In the subsequent section[3l based 
on this ansatz, the backreaction onto the motion of the fluid using the full 
current will be calculated, and it is shown that this yields a different result 
than that obtained by the effective-action method. Afterwards, the failure of 
the effective-action technique is discussed in more detail in section |4| The 
cutoff dependence of the pseudo-energy-momentum tensor ([9]) is addressed in 
section [5] As a simple example, we consider the influence of the backreaction 
contribution on a static quasi-lD condensate in section [Hi 



2 Excitations in Bose-Einstein Condensates 
2.1 Particle-number-conserving mean-field expansion 

In the s-wave scattering approximation, a dilute many-particle system of inter- 
acting bosons is described, on a "microscopic" level corresponding to distance 
scales much larger than the true range of the interaction potential, by the 
field operator equation of motion in the Heisenberg picture (we set from now 
on ?i = m = 1) 

^i-A^ + (11) 
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In the limit of many particles iV 3> 1, in a finite trap at zero temperature with 
almost complete condensation, the full field operator ^ can be represented in 
terms of the particle-number-conscrving mean-field ansatz |301 131] 

#= (V^e+X + C)i^"'/'. (12) 

Here, the order parameter ■i/'c ~ 0{\/N) [note that ipc ^ (^) , as opposed to the 
•f' in Eq. The one-particle excitations arc denoted x = 0{N^), where one- 
particle here means that the Fourier components of x a-rc linear superpositions 
of annihilation and creation operators of quasiparticles ak and a\,, cf. Eq. (jSSp 
below. The remaining higher-order, multi-particle corrections arc described by 
C, = 0{1/^/N). The above mean- field ansatz can be derived in the dilute-gas 
hmit by formally setting g = 0{l/N) [HIISIISIISI]; we shall use this formal 
definition of the dilute-gas limit in what follows. The dilute-gas limit should 
be compared and contrasted with the usual thermodynamic limit, in which 
the density and particle interaction remains constant, while the size of the 
(trapped) system increases with N oo, adjusting the harmonic trapping 
potential Vtrap correspondingly (in D spatial dimensions, the thermodynamic 
limit corresponds to keeping Nui^ constant for N oo, where w is the 
geometric mean of the trapping frequencies ^ ) . In the presently used dilute- 
gas limit, on the other hand, the trapping potential remains constant, but the 
interaction and the density change. The advantage of the limit gN constant is 
that in this limit we have a well-defined prescription to implement the mean- 
field approximation, keeping one power of g for each factor of TV, cf. |36| . 



d 



2.2 Gross-Pitaevskii and Bogoliubov-de Gennes equations 

Insertion of Eq. (|12p into pT|) yields to 0{N) the Gross-Pitaevskii equation 
|35| for the order parameter ip^ 

V'c= (-^V2 + yt,,p+,g|^.eP + 2g(x^x))^c + ff(x')V^c*- (13) 

The Bogoliubov-de Gennes equations [37| for the one-particle fluctuations x 
are obtained to 0(iV°) 

'Wt^ = (-^ + Ftrap + 2.g|7^ep) X + Q^X^ , (14) 

Finally, the time evolution of the remaining higher-order corrections in the 
expansion ([72]) . C, = 0{\/^fN), neglecting the 0{1/N) terms, is given by: 



+2.9(X^X - (xtx))^c + 5(x' - {X^))rc ■ (15) 



The Gross-Pitaevskii equation in the form ([13]) ought to be compared with 
the simple-minded form of The additional terms 2g{x^x) ^J^d g{x^) in the 
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Gross-Pitacvskii equation in the form of Eq. above ensure that the expec- 
tation value of the muhi-particle operator, ( = 0{1/\/N), vanishes in leading 
order, {0 = 0{1/N). Without these additional terms, the mean-field expan- 
sion (ini) would still be vahd with C = 0(1///V), but without (C) = 0{l/N). 
The proper incorporation of the so-called "anomalous" fluctuation average 
(x^) (the "normal" fluctuation average is (x^x)) into the description of Bose- 
Einstein condensates has also been discussed from various points of view in 



3 Quantum Backreaction 

3.1 Calculation of backreaction force from microscopic physics 

The observation that the Gross-Pitaevskii equation yields an equation 
correct to leading order 0{\/N), using either I'lpd'^ or iV'cP + 2 (x'x) the 
first line of (fT3|) . hints at the fact that quantum backreaction effects corre- 
spond to next-to-leading, i.e. quadratic order terms in the fluctuations and 
cannot be derived ab initio in the above manner without additional assump- 
tions. Therefore, we shall employ an alternative method: In terms of the exact 
density and current given by 

^, = , j = 1 vif- - H.c.\ , (16) 



2i 

the time-evolution is governed by the equation of continuity for g and an 
Euler type equation for the current j. After insertion of Eq. pTjl . we find that 
the equation of continuity is not modified by the quantum fluctuations but 
satisfied exactly (i.e., to all orders in or h) 

^f?+V-j = 0, (17) 

in accordance with the U{1) invariance of the Hamiltonian and the Noether 
theorem, cf. [28| . However, if we insert the mean-field expansion (|12p and write 
the full density as a sum of condensed and non-condensed parts 

Q = gc + {x^x)+o{i/VN), (18) 

with Qc — iV'cP, we find that neither part is conserved separately in general. 
Note that this split requires {() = 0{1/N), i.e., the modifications to the Gross- 
Pitaevskii equation (jl3p discussed above. Similarly, we may split up the full 
current [with QcVc = 3(V'c^'0c)] 

3 = QcVc + ^ (x^Vx - H.c.) + 0{1/Vn) , (19) 

and introduce an average velocity v via j ~ gv. This enables us to unambigu- 
ously define the quantum backreaction Q as the following additional contri- 
bution in an equation of motion for j analogous to the Euler equation: 
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= /ci(j, 9) + Q + 0(1/ Vn) , (20) 
where the classical force density term 

/ciW, 9) = -v[V- (gv)] - g{v ■ V)v + gV (^^^^ - ^trap - gg^ .(21) 

Here, "classical" means that the force density contains no explicit quantum 
fluctuation terms (i.e., only those absorbed in the full density and full current), 
and in addition just the "quantum pressure", which already occurs on the 
mean-field level. Formulation in terms of the conventional Euler equation, 
i.e., using a convective derivative of the velocity defined by j = pv giving the 
acceleration, yields 

I + . ^ V) . . - V + - 1 ^) + 5 + o(N-"") , (22) 

The quantum backreaction force density Q can now be calculated by com- 
paring the two equations above and expressing dj /dt in terms of the field 
operators via Eqs. (fTTI) and (fTB|l 

-T-J = - /•f't - (V^*'t)V# + H.c. 

- vVtrap - l^gH^^ f^'-) . (23) 

After insertion of the mean-field expansion (|12[) , we obtain the leading con- 
tributions in the Thomas-Fermi limit 

Q = V-{y®j-^+j-^®v - g^v (g) v) 

-^V.((Vx^)®Vx + H.c.) , (24) 

with g-^ = (x^x) and j'^ = 9 (x^'^x)- Under the assumption that the relevant 
length scales A for variations of, e.g., g and g, are much larger than the healing 
length ^ = (.gg)~^/^, we have neglected terms containing quantum pressure 
contributions V^£> and [V^)]^, which amounts to the Thomas- Fermi or local- 
density approximation. These contributions would, in particular, spoil the 
effective (local) geometry in Eq. ([S]) (the inclusion of the quantum pressure to 
derive a "nonlocal metric" has been discussed in [40]). 



3.2 Comparison with effective-action technique 

To compare the expression for the backreaction force density derived from 
the full dynamics of ([24]) with the force obtained from Eq. (|10p . we have 
to identify the scalar field (p and the contravariant metric g^''. We already 
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know that phonon modes with wavelength A 3> ^ are described by the action 
in Eq. ^ in terms of the phase fluctuations (j) provided that g'"^ is given by 
([31), where A^. = Cg/g and y^—g = c^/g'^ = p^/cs in three spatial dimen- 
sions. The density fluctuations 5g are related to the phase fluctuations cf) via 
Sg= -g-^{d/dt + v ■V)cj). 

The variables 77 = {pb,</'b} or alternatively 77 = {pb, V(/)b} = {pbji'b} in 
(dni) are then defined by the expectation values of density and phase operators 
according to 

<P = (<^) + = 0b + , 

Q = {q) + ^ gh + Sg. (25) 

The phase operator can formally be introduced via the following ansatz for 
the full field operator 

# = e'*V^. (26) 

Since ^ and g do not commute, other forms such as 'F ~ \fge^^ would not 
generate a self-adjoint ^ (and simultaneously satisfy 'i''^'^ = g). Note that, 
in contrast to the full density which is a well-defined and measurable quan- 
tity, the velocity potential, ^, is not [IT]. This can be seen as follows. The 
commutator between density and phase operators 

[g{r)Mr')]^i5[r-r'), (27) 

yields, if one takes, on both sides, matrix elements in the number basis of the 
space integral over r for a given volume V in which r' is situated, 

[N ~ N'){N\^{r')\N') = zW, (28) 

where N and N' are two possible values for the number of particles in V. This 
is an inconsistent relation for N, N' positive semidefinite and discrete (which 
the very existence of particles requires), most obviously for N = N' . The 
commutator therefore makes sense only if it is understood to be effectively 
coarse-grained over a sub-volume V with large enough number of particles 
N ^ 1, such that the inconsistency inherent in (|28p becomes asymptotically 
irrelevant. It cannot be defined consistently locally, i.e., in arbitrarily small 
volumina, where there is just one particle or even none, or when the number 
fluctuations in larger volumina V are large, such that the probability to have 
a very small number of particles in V is not negligible [42j . 

The action in terms of the total density g and the variable ^ reads (ne- 
glecting the quantum pressure term, i.e., in the Thomas-Fermi limit) 

C^-g (^<^+ ^C^<^)') - ^[^'l - Vtr.„g, (29) 

with e[g] denoting the internal energy density. The quantum corrections to the 
Bernoulli equation up to second order in the fiuctuations, using the effective- 
action method in Eq. (fTO)) are then incorporated by writing (the background 
density here equals the full density, g ~ gh) 
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^•^b + \{^<l^^.f + h[Q] - ^ = , (30) 
where h[Q\ ~ de/dg + Vtrap- We obtain, using ([5]), 

Clearly, taking the gradient of this result we obtain a backreaction force which 
markedly differs from the expression (|24|) derived in the previous subsec- 
tion. Moreover, it turns out that the backreaction force density in Eq. 
contains contributions which are not part of the expectation value of the 
pseudo-energy-momentum tensor, {T^i,). For example, the phonon density 
Q-^ contains (((5^)^)ren (where (. ..)ren means that the divergent c-number 
Xx^ — x^x — '5(0) has been subtracted already) which is part of (T^i/)ren, 
but also ((/)^)ron which is not. [Note that (0^)ron cannot be cancelled by the 
other contributions.] The expression ((Vx^) ® Vx + H.c.) in the last line 
of Eq. p4|) contains (V(/) (g) V(/))ren which does occur in {Tfj,^)i-en, but also 
{^Sg (g) V(5^)ron, which does not. One could argue that the latter term ought 
to be neglected in the Thomas-Fermi or local-density approximation since it 
is on the same footing as the quantum pressure contributions containing V^p 
and [Vg]'^ [which have been neglected in ([M)) ]. but it turns out that this 
expectation value yields cutoff dependent contributions of the same order of 
magnitude as the other terms, see section [5] below. 



4 Failure of Effective- Action Technique 

After having demonstrated the failure of the effective-action method for de- 
ducing the quantum backreaction, let us study the reasons for this failure in 
more detail. The full action governing the dynamics of the fundamental field 
'F reads 

= i<p*^q^ (^-iv' + Krap + l'^'*^') f'. (32) 

Linearization according to = i/^c -I- X yields the effective second-order action 
generating the Bogoliubov-de Gennes equations ^4]) 

'C^ff = *X*^X - X* (-^V^ + Ftrap + 2g\i^l\^ x-\ [5(^c*)'x' + H.c] . 

(33) 

If we start with the action ([29|) in terms of the total density g and the nonfun- 
damental variable the quantum corrections to the equation of continuity 
(5 / 50b are reproduced correctly but the derived quantum backreaction con- 
tribution to the Bernoulli equation, SAcsl ^Q\i^ and therefore the backreaction 
force, is wrong. 
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One now is led to the question why the effective-action method works for 
the fundamental field 'Jf and gives the correct expression for the backreaction 
force, but fails for the non-fundamental variable The quantized fundamen- 
tal field ^ satisfies the equation of motion pT|) as derived from the above 
action and possesses a well-defined linearization via the mean-field expansion 
(fT2|) . One of the main assumptions of the effective-action method is a similar 
procedure for the variable ^, i.e., the existence of a well-defined and lineariz- 
able full quantum operator <1> satisfying the quantum Bernoulli equation (for 
large length scales), cf. Eq. ([50]) 

y>+\{v-l>f + h[e]lQ. (34) 

The problem is that the commutator of g and ^ at the same position diverges, 
cf. Eq. (P7|) . and hence the quantum Madelung ansatz in Eq. ([25]) is singular. 
As a result, the above quantum Bernoulli equation is not well-defined (in con- 
trast to the equation of continuity), i.e., insertion of the quantum Madelung 
ansatz in Eq. ([26]) into Eq. (fTT|) generates divergences [41] . 

In order to study these divergences by means of a simple example, let us 
consider a generalized Bose-Hubbard Hamiltonian [IS] , which considers bosons 
sitting on a lattice with sites i, which can hop between nearest neighbor sites 
and interact if at the same site: 

^ = - f E ^^l^^ + H-c.) + E [l^^^^ + i^') ' (35) 

<ij> i 

where is the annihilation operator for bosons at a given lattice site i, and 
< ij > denotes summation over nearest neighbors; rij = '^j'l^i is the so-called 
filling factor (operator), equal to the number of bosons at the lattice site 
i. The quantities Pi multiplying the filling factor depend on the site index. 
In the continuum limit, the lattice Hamiltonian (|35p generates a version of 
Eq. (|TT]) . Setting a^^^i'{xi) = where the 'i'{xi) are the continuum field 
operators and a is the lattice spacing taken to zero (we consider for simplicity 
a simple cubic lattice in D spatial dimensions), the effective mass is given by 
1/m* = aa^: The bosons moving through the lattice obviously become the 
heavier the smaller the hopping amplitude a becomes at given a. The coupling 
constant is determined by g = 70^, and the trap potential is governed by 
Vtre,pix)=g/2 + (3{x)~a. 

On the other hand, inserting the quantum Madelung ansatz employing a 
phase operator, Eq. (j26p in its lattice version, the problem of operator ordering 
arises and the (for the Bernoulli equation) relevant kinetic energy term reads 

H^ = \Y. VM^^^TT) (V<^)2 + H.c. , (36) 

i 

with the replacement fii + I instead of fii being one effect of the non-commu- 
tativity. In the superfluid phase with large filling n 3> 1, we therefore obtain 
the following leading correction to the equation of motion 
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at 2 n 16 n J 

The BernouUi equation in its quantum version thus receives corrections de- 
pending on microscopic details Hkc the fiUing of a particular site and is not of 
the (conjectured) form 

By means of this simple example, we already see that the various limiting 
procedures such as the quantization and subsequent mean-field expansion, 
the variable transformation {^* ^W) ^ (f?, ^^), and the linearization for small 
fluctuations, as well as continuum limit do not commute in general - which 
explains the failure of the effective-action method for deducing the quantum 
backreaction. The variable transformation (<P*,'F) ^ (g,^) is applicable to 
the zeroth-order equations of motion for the classical background as well as 
to the first-order dynamics of the linearized fluctuations - but the quantum 
backreaction is a second-order effect, where the aforementioned difficulties, 
such as the question of the choice of fundamental variables and their operator 
ordering, arise. 



5 Cutoff Dependence of Effective Action 

Another critical issue for the applicability of the effective-action method is 
the UV divergence of {T^^). Extrapolating the low-energy effective action in 
Eq. ^ to large momenta k, the expectation values {Sg^) and {(f)^) entering 
would diverge. For Bose-Einstein condensates, we may infer the deviations 
from Eq. ([5|) at large k from the Bogoliubov-de Gennes equations (fT4ll . Assum- 
ing a static and homogeneous background (which should be a good approxima- 
tion at large k), a normal-mode expansion yields a Bogoliubov transformation 
between the bare bosonic operators Xk and the quasiparticle operators a^, a\^: 



Xk 




(38) 



where the frequency ojk is determined by the Bogoliubov dispersion relation 
for the dilute Bose gas, = ggk^ + k^/A. The above form of the Bogoliubov 
transformation results, after inversion, in the usual phonon quasiparticle op- 
erators at low momenta, and gives Xfc = Ofc at fc — > oo, i.e., the quasiparticles 
and the bare bosons become, as required, identical at large momenta. 

Using a linear dispersion oc k^ instead of the full Bogoliubov dispersion, 
expectation values such as (x'x) would be UV divergent, but the correct 
dispersion relation implies Xk ~ o,k9g/k^ + ^fc for large k^, and hence (x^x) 
is UV finite in three and lower spatial dimensions. Thus the healing length ^ 
acts as an effective UV cutoff, = 1/^. 

Unfortunately, the quadratic decrease for large k in Eq. psp . giving asymp- 
totically Xk ^ o,k9g/f^^ + flfei is not sufficient for rendering the other expec- 
tation values (i.e., apart from g^ and j ) in Eq. ([M]) UV finite in three spatial 
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dimensions. This UV divergence indicates a failure of the s-wavc pseudo- 
potential gS^{r — r') in Eq. (fTTI) at large wavenumbers k and can be elim- 
inated by replacing gS^{r — r') by a more appropriate two-particle interac- 
tion potential Vint{r — r'), see [35] ■ Introducing another UV cutoff wavenum- 
ber fc^"' related to the breakdown of the s-wave pseudo-potential, we obtain 

((Vx^) «) Vx + H.C.) ~ ;,cut g^^^ ^^2^ _ f^cut^ 

In summary, there are two different cutoff wavenumbers: The first one, 
/c|"', is associated to the breakdown of the effective Lorentz invariance (change 
of dispersion relation from linear to quadratic) and renders some - but not all 
- of the naively divergent expectation values finite. The second wavenumber, 
/Cg"*, describes the cutoff for all (remaining) UV divergences. In dilute Bose- 
Einstein condensates, these two scales are vastly different by definition. Be- 
cause the system is dilute, the inverse range of the true potential fc™* = l/rp, 
must be much larger than the inverse healing length. Thus the following con- 
dition of scale separation must hold: 

l,CUt _ I,CUt ^ 7,CUt _ LCUt /qQ\ 

KUV — ^ 'S^ — ft^Lorontz • 

In terms of the length scales governing the system, using that fc|^* = y/iirasP, 
in the dilute gas we must have the following condition fulfilled, inagrQ <^ d'^, 
where d = g~^^^ is the interparticle separation. Note that the opposite scale 
separation relation, kf^^^^^^ ^ fc^y is very unnatural in the sense that every 
quantum field theory which has the usual properties such as locality and 
Lorentz invariance must have UV divergences (e.g., in the two-point function). 

The renormalization of the cutoff-dependent terms is different for the two 
cases: The fcg^^-contributions can be absorbed by a ^(-independent renormal- 
ization of the coupling g [34l [36] , whereas the fc™*-contributions depend on 
the density in a nontrivial way and thus lead to a quantum renormalization of 
the effective equation of state. We supply an example for this renormalization 
in the section to follow. 



6 Static Example for the Backreaction Force 

In order to provide an explicit example for the quantum backreaction term 
in Eq. ((24]) , without facing the above discussed UV problem, let us consider a 
quasi-one-dimensional (quasi-lD) condensate |44l|45], where all the involved 
quantities are UV finite. In a quasi-lD condensate the perpendicular harmonic 
trapping uj±_ is much larger than the axial trapping lj, such that the conden- 
sate assumes the shape of a strongly elongated cigar with all atomic motion in 
the perpendicular direction frozen out, ujj_ being much larger than the mean 
energy per particle. 

In accordance with general considerations [46], the phonon density is 
infrared (IR) divergent in one spatial dimension, therefore inducing finite-size 
effects, i.e., a dependence of the various quantities of interest on the system 
size. Nevertheless, in certain situations, we are able to derive a closed local 
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expression for the quantum backrcaction term Q: Let us assume a completely 
static condensate t; = in effectively one spatial dimension, still allowing for a 
spatially varying density g and possibly also coupling g. Furthermore, since we 
require that spatial variations of g and g occur on length scales A much larger 
than the healing length (Thomas- Fermi approximation), we keep only the 
leading terms in ^/A <C 1, i.e., the variations of g and g will be neglected in the 
calculation of the expectation values. In this case, the quantum backrcaction 
term Q simplifies considerably and yields (in effectively one spatial dimension, 
where g = giD and g = giD now both refer to the ID quantities) 

Q = -V ((Vxt)Vx) - {g'-g {2x^x + ix^f + x')) 



-V 



o(,e/\') . (40) 



It turns out that the IR divergences of 2(x^x) a-nd ((x^)^ + x^) cancel each 
other such that the resulting expression is not only UV but also IR finite. 
Note that the sign of Q is positive and hence opposite to the contribution of 
the pure phonon density (x'x), which again illustrates the importance of the 
"anomalous" term ((x^)^ +X^)- 

A possible experimental signature of the quantum backreaction term Q 
calculated above, is the change incurred on the static Thomas-Fermi solution 
of the Euler equation for the density distribution (cf. [551 136j ) 



= ^i^i^ + VKI^ + oil/VJ^) , (41) 

with /i denoting the (constant) chemical potential. The classical [0{N)] den- 
sity profile gd = (/U — Vtrap)/5iD acquires nontrivial quantum [©(A'^^)] cor- 
rections gq = -^/i — Vtrap/27r, where the small parameter is the ratio of the 
interparticle distance 1/ g = 0{1/N) over the healing length ^ = 0{N^). Note 
that the quantum backreaction term gg in the above split g = gd + gq should 
neither be confused with the phonon density g^in g = gc + g^ (remember that 
g^ is IR divergent and hence contains finite-size effects) nor with the quantum 
pressure contribution cx V^y^ in the Euler type Eq. ([22]) . 

Evaluating the change AR of the Thomas-Fermi size (half the full length), 
where ^ = Vtrap, of a quasi-lD Bosc-Einstein condensate induced by back- 
reaction, from Eq. (|¥T|) we get Z\i? = — 2^^/^((jjj^/(jj2)as. Here, the quasi- 
lD coupling constant ^id is related to the 3D s-wave scattering length 
Us and the perpendicular harmonic trapping uj± by ^id = 2asU!± (provided 
<C = l/y/uj± [44]). In units of the classical size Rd ~ {3asNuj±/Ljfy/^ , 
we have 



AR _ 1 / 1 y f^^± a. 



1/3 / ^ 2/3 

(42) 
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where ~ \/ ^JuJl describes the longitudinal harmonic trapping length. In 
quasi- ID condensates, backreaction thus leads to a shrinking of the cloud rel- 
ative to the classical expectation - whereas in three spatial dimensions wc 
have the opposite effect [H I36j . In one dimension, we thus obtain a soften- 
ing of the quantum renormalized equation of state of the gas; conversely, in 
three spatial dimensions the effective equation of state becomes stiffer due to 
quantum fluctuations. 

For reasonably realistic experimental parameters, the effect of quantum 
backreaction on the equation of state should be measurable; for N ~ 10^, 
ui_l/uJz ~ 10^, and a^/a^ ~ 10~^, we obtain \AR/Rc\\ ~ 1%. 



7 Conclusion 

By explicit analysis of the analytically tractable case of a dilute Bose gas 
in the mean-field approximation, we have demonstrated the following. Even 
given that the explicit form of the quantum backreaction terms depends on 
the definition of the classical background, the effective-action method does not 
yield the correct result in the general case (which is, in particular, independent 
of the choice of variables). The knowledge of the classical (macroscopic) equa- 
tion of motion - such as the Bernoulli equation - may be sufficient for deriving 
the first-order dynamics of the linearized quantum fiuctuations (phonons) , but 
the quantum backreaction as a second-order effect cannot be obtained with- 
out further knowledge of the microscopic structure (which rcfiects itself, for 
example, in the operator ordering imposed). It is tempting to compare these 
findings to gravity, where wc also know the classical equations of motion only 

which - in analogy to the Bernoulli equation - might yield the correct first- 
order equations of motion for the linearized gravitons, but perhaps not their 
(second-order) quantum backreaction. Another potentially interesting point of 
comparison is the existence of two different high-energy scales - one associated 
to the breakdown of Lorentz invariance fc™* = kf^^^^^^ and the other one, 
^uv = '^^s"': to the UV cutoff introduced by the true interaction potential 
range. The question then poses itself whether one of the two cutoff scales (or 
in some sense both of them) correspond to the Planck scale in gravity. 

The dominant ©(^/A) quantum backreaction contributions like those in 
Eq. (|¥T|) depend on the healing length as the lower UV cutoff and hence can- 
not be derived from the low-energy effective action in Eq. ([5]) using a covari- 
ant (i.e., cutoff independent) regularization scheme, which does not take into 
account details of microscopic physics (represented, for example, in the quasi- 
particle dispersion relation). Note that the leading 0(^/A) quantum correction 
to the pressure could be identified with. &. cosmologica-1 term, (X'^i/) — yl g/n^ 
in Eq. ([TD|), provided that the cosmological "constant" A is not constant but 
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depends on g and g. Note that in general relativity, the Einstein equations 
demand A to be constant, due to the equivalence principle and the resulting 
requirement that the metric be parallel-transported, V^g^^ = 0. 

As became evident, the knowledge of the expectation value of the pseudo- 
energy- momentum tensor {T^^) is not sufficient for determining the quantum 
backreaction effects in general. Even though {T^^) is a useful concept for de- 
scribing the phonon kinematics (at low energies), we have seen that it does not 
represent the full dynamics of the fluid dynamical variables defined in terms of 
the fundamental quantum field ^ . Related limitations of the classical pseudo- 
energy-momentum tensor, in particular the background choice dependence of 
the description of the second-order effect of the exchange of energy and mo- 
mentum between excitations and that background, and the resulting form of 
the conservation laws, have been discussed in |27| . 

In general, the quantum backreaction corrections to the Euler equation 
in Eq. ([22]) cannot be represented as the gradient of some local potential, 
cf. Eq. (|24p . Hence they may effectively generate vorticity and might serve as 
the seeds for vortex nucleation from the vortex vacuum. 

In contrast to the three-dimensional case (see, e.g., [35l [36]), the quan- 
tum backreaction corrections given by Eq. (|40p diminish the pressure in con- 
densates that can be described by Eq. pT|) in one spatial dimension (quasi- 
ID case). This is a direct consequence of the so-called "anomalous" term 
{(x^Y + X'^) in Eq. PO)) . which - together with the cancellation of the IR di- 
vergence - clearly demonstrates that it cannot be neglected in general. We 
emphasize that even though Eqs. (|40)) - (|42p describe the static quantum back- 
reaction corrections to the ground state, which can be calculated by an alter- 
native method [36] as well, the expression in Eq. ([24]) is valid for more general 
dynamical situations, such as rapidly expanding condensates. Quantum back- 
reaction can thus generally not be incorporated by rewriting the Euler equa- 
tion in terms of a renormalized chemical potential. While the static quantum 
backreaction corrections to the ground state can be absorbed by a redefi- 
nition of the chemical potential //(p) determining a quantum renormalized 
(barotropic) equation of state p{q)^ this is not possible for the other terms in 
Eq. (|24p . like the quantum friction-type terms depending on j ® v. 

We have derived, from the microscopic physics of dilute Bose-Einstein 
condensates, the backreaction of quantum fluctuations onto the motion of the 
full fluid and found a quantum backreaction force that is potentially exper- 
imentally observable in existing condensates. We observed a failure of the 
effective- action technique to fully describe the backreaction force in Eq. ([24|), 
and a cutoff dependence of backreaction due to the breakdown of covariance 
at high energies. Whether similar problems, in particular the question of the 
correct choice of the fundamental variables and the related operator ordering 
issues, beset the formulation of a theory of "real" (quantum) gravity remains 
an interesting open question. 
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